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Abstract 



In this paper, we study the unconditional convergence and error estimates 
of a Galerkin-mixed FEM with the linearized semi-implicit Euler time-discrete 
scheme for the equations of incompressible miscible flow in porous media. We 
prove that the optimal L'^ error estimates hold without any time-step (conver- 
gence) condition, while all previous works require certain time-step condition. 
Our theoretical results provide a new understanding on commonly-used lin- 
earized schemes. The proof is based on a splitting of the error function into 
two parts: the error from the time discretization of the PDEs and the error from 
the finite element discretization of corresponding time-discrete PDEs. The ap- 
proach used in this paper is applicable for more general nonlinear parabolic 
systems and many other linearized (semi)-implicit time discretizations. 

1 Introduction 

Incompressible miscible flow in porous media is governed in general by the following 
system of equations: 



where p is the pressure of the fluid mixture, u is the velocity and c is the concen- 
tration of the first component; k{x) is the permeability of the medium, /x(c) is the 
concentration-dependent viscosity, <1> is the porosity of the medium, and q^ are 

'Department of Mathematics, City University of Hong Kong, Kowloon, Hong Kong. The work 
of the authors was supported in part by a grant from the Research Grants Council of the Hong 
Kong Special Administrative Region, China (Project No. CityU 102005) 

Email address: buyangli2@student.cityu.edu.hk, maweiw@math.cityu.edu.hk 







(1.3) 



given injection and production sources, c is the concentration of the first component 
in the injection source, and D{\i) = [Dij{u)]iixd is the diffusion-dispersion tensor 
which may be given in different forms (see [H [5] for details) . We assume that the 
system is defined in a bounded smooth domain Q in R'^ {d = 2,3), for t £ [0,T], 
coupled with the initial and boundary conditions: 

u-n = 0, D(u)Vc-n = for x G t e \0,T], 

(1.4) 

c(x,0)=co(x) for X £ Q,. 

The above system was investigated extensively in the last several decades [3 El 
[T7 1 ll8 1 l20] due to its wide applications in various engineering areas, such as reservoir 
simulations and exploration of underground water, oil and gas. Existence of weak 
solutions of the system was obtained by Feng |20) for the 2D model and by Chen 
and Ewing [6] for the 3D problem. Existence of semi-classical/classical solutions 
is unknown so far. Numerical simulations have been done extensively with various 
applications, see [10^ IT2| [29] and the references therein. Optimal error estimates of 
a standard Galerkin-Galerkin method for the system ()l.ip - (jl.4p in two-dimensional 
space was obtained first by Ewing and Wheeler [19] roughly under the time-step 
condition r = o{h). In their method, a linearized semi-implicit Euler scheme was 
used in the time direction and Galerkin FEM approximation was used both for the 
concentration and the pressure. Later, a Galerkin-mixed finite element method was 
proposed by Douglas et al |11] for this system, where a Galerkin approximation was 
applied for the concentration equation and a mixed approximation in the Raviart- 
Thomas finite element space was used proposed for the pressure equation. A lin- 
earized semi-implicit Euler scheme, the same as that used in [19], was applied for the 
time discretization. Optimal error estimates were obtained under a similar time-step 
condition r = o(/i). In [28j, the authors studied the system in two-dimensional space 
by the linearized backward Euler scheme and a streamline-upwind-Petrov-Galerkin 
method combined with a post-process technique on the velocity. Error estimates 
with quasi-optimal rates were derived. There are many other numerical methods in 
literature for solving the equations of incompressible miscible flow, such as |35[ [38] 
with an ELLAM in two-dimensional space, ^36j with an MMOC-MFEM approxi- 
mation for the 2D problem, |33[ [T3] with a characteristic-mixed method in two and 
three dimensional spaces, respectively, and [26^127] with a collocation- mixed method 
and a characteristic-collocation method, respectively. In all those works, error es- 
timates were established under certain time-step conditions. Moreover, linearized 
semi-implicit schemes have also been analyzed for many other nonlinear parabolic- 
type systems, among which many require certain time-step restrictions when strong 
nonlinearities were involved, see [II |21 13 E [El [13 El [^ |2ll |25l [371 [38] and the 
references therein. Such time-step restrictions may result in the use of a very small 
time step and extremely time-consuming in practical computations. The problem 
becomes more serious when a non-uniform mesh is used. However, we believe that 
those time-step restrictions are not necessary in most cases. 

In this paper, we analyze the linearized semi-implicit Euler scheme with a pop- 
ular Galerkin-mixed finite element approximation in the spatial direction for the 
system (|l.ip - (|1.4p . We establish optimal error estimates almost without any 
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time-step restriction (or when h and r are smaller than some positive constants, 
respectively). This provides a new understanding on the commonly-used linearized 
schemes and clears up the misgivings for the time-step size in practical computa- 
tions. Our theoretical analysis is based on an error splitting proposed in [23] for a 
Joule heating system with a standard Galerkin FEM. By introducing a correspond- 
ing time-discrete system, we split the numerical error into two parts, the error in 
the spatial direction and the error in the temporal direction. The rigorous analysis 
of the regularity of the solution to the time-discrete PDEs is a key to our approach 
because of the strong coupling of the system. With such proved regularity, we obtain 
optimal and r-independent L^-error estimates of the Galerkin-mixed FEM for the 
time-discrete PDEs. 

The rest of the paper is organized as follows. In Section 2, we introduce the 
linearized semi-implicit Euler scheme with a Galerkin-mixed approximation in the 
spatial direction for the system ()l.ip - (ll.4p and present our main results. In Section 
3, we present a priori estimates of the solution to the corresponding time-discrete 
system and the error estimate of the linearized scheme. Optimal error estimates of 
the fully discrete scheme in the norm are given in Section 4. Conclusions are 
drawn in the last section. Since numerical simulations on the underlying system 
have been done extensively, we will not provide any numerical results in this paper. 



2 The Galerkin-mixed FEM and the main results 

For any integer m > and 1 < p < oo, let VK™'^ be the Sobolev space of functions 
defined on equipped with the norm 



i/3|<m 

where 



dm 



1 



5xf • • • dx'^/ 

for the multi-index /3 = (/3i, • • • ,Pd), Pi > 0, • • • , /3d > 0, and = /3i H h Pd- 

For any integer m > and < a < 1, let C"""*"" denote the usual Holder space 
equipped with the norm 

\f}\<m \p\=k 

Let / = (0,T). For any Banach space X, we consider functions g : I ^ X and 
define the norm 



\\9\\lp{I;X) — ' 



T \ - 

^ \\git)\f^dty , l<p<oo. 



esssuptg^||5f^tj||x, p = oo 
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In addition, we define Lq as the subspace of Lp consisting of functions in Lp whose 
integral over Q are zeros. Finally, we denote by H the space of vector-valued func- 
tions f £ L"^ X L'^ X L'^ such that V ■ f £ L'^. 

Let TTh be a regular subdivision of fl into triangles Tj, j = 1, • • • ,M, in (or 
tetrahedra in M^), with fi/j = UjTj and denote hy h = maxi<j<jv/{diam Tj} the 
mesh size. For a triangle Tj at the boundary, we define Tj as the extension of Tj to 
the triangle with one curved edge (or a tetrahedra with one curved face in M^). For 
a given subdivision of il, we define the finite element spaces: 

Sh = {wh G L\n) : Wh\f. is linear for each element Tj £ tth}, 
= {wh £ C^{^) : Wh\f. is linear for each element Tj £ TTh}- 

Let Hh be the subspace of H, as introduced by Raviart and Thomas |3T1 [M] such 
that divi/j £ Sh for £ Hh- 

In the rest part of this paper, we assume that the solution to the initial-boundary 
value problem (ll.lj) - (ll.4p exists and satisfies 

+ ||u||icx)(/.//2) + ||ut||£^2(/.tyl,3/2) + ||c||ioo(/.(4/2,s) 
+ l|ct||L°o(7;L«) + l|ct||L''(/;VKi.'') + l|Ctt||L^(/;H/-l>'=) < C (2.1) 

for some s > d and 

ll/lbi, II/IIhi <C. (2.2) 
Correspondingly, we assume that the permeability k £ C'^{Q) and satisfies 

< k{x) < ko for x £ n, (2.3) 
the concentration-dependent viscosity ^ £ C^(M) and satisfies 

Mo ^ ^ /^(s) ^ /^o for s G R, (2.4) 

for some positive constant fiQ. Moreover, we assume that the diffusion-dispersion 
tensor is given by D{u) = ^dmI+D*{u), where dm > 0, -D*(u) = di{u)I +d2{u){u0 
u) is symmetric and positive definite and d^.^.D £ L°° [3]. For the initial-boundary 
value problem (|l.ip - (|1.4p to be well-posed, we require 

/ /dx= / g^dx. (2.5) 

Jo, 

Let {tn]n=o ^ uniform partition of the time interval [0, T] with r = T /N and 
denote 

p"-=p{x,tn), u" = u(x,tn), = c(x, t„) . 

For any sequence of functions {/"}^=0' define 

fn+l £71 

^ J J_ 

T 
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The fully discrete mixed finite element scheme is to find Pf^ E 5'/j/{constant}, 
Uf^ S Hfi and G Vh, n = 0,1, ■ ■ ■ , N, such that for all Vh G Hh, (fh £ Sh and 

— t^r'> ^/.V -fn"^'' V • (2-6) 



(V • C/^"+\ V^,) = (g^-g^, v^/.), (2.7) 

($ACr\ <t>h) + {d{U]:+^)VCI+\ V0,) 

+ (c/,"+i • VC+i, 0,) = (cq' - Cl+\P, 0,) (2.8) 

where the initial data is the Lagrangian piecewise linear interpolation of . 

In this paper, we denote by C a generic positive constant and by e a generic 
small positive constant, which are independent of n, h and r. We present our main 
results in the following theorem. 

Theorem 2.1 Suppose that the initial-boundary value problem M.1]) - ^^L^ has a 

unique solution (p, u, c) which satisfies Then there exist positive constants ho 

and To such that when h < ho and r < tq, the finite element system i2.6\) - [278\) 
admits a unique solution {PJ^, UJ^,Cf^), n = 1, • • • , A^, which satisfies that 

max IIP," + max ||^," - u"||^2 + max ||C " c^l^ < C{t + h\ 

l<n<N l<n<N l<n<N 

(2.9) 

We will present the proof of Theorem 2.1 in the next two sections. The key to 
our proof is the following error splitting 

llTTTi ,,n|| ^ ll„i|| I llrr" TTi^ll 

||ty, - u 11^2 < ||e„||i2 + ||(y -tV, ||l2, 

1 1 P" ^"-11 „ ^ 1 1 1 1 „ _L 1 1 P" 1 1 „ 

- c 11^2 < lle^ ||l2 + ||C -L^Wui, 

where 



e" 


pn 






= [/" 


-u" 


e" 




-c". 



and (-P", U"',C^) is the solution of a time-discrete system defined in next section. 

3 Error estimates for time-discrete system 

We define the time-discrete solution (P", C/"',C") by the following elliptic system: 

^n+l ^ __^VP"+\ (3.1) 
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V • =q^ - qP, (3.2) 

$AC"+^ - V • (D(;7"+i)VC"+^) + [/"+^ • VC"+^ = cq^ - C"+ig^, (3.3) 

for X G $7 and t G [0, T], with the initial and boundary conditions 

[/"+!. n = 0, i:>([/"+i)VC"+i • n = for x G 917, tG[0,r], 
C''(x) = Co(x) for X G il, 

In this section, we prove the existence and uniqueness and suitable regularity of 
the solution of the above time-discrete system. 

Theorem 3.1 Suppose that the initial-boundary value problem hl.l\) - ^L^ has a 
unique solution {p, u, c) which satisfies ^2.1\) . Then there exists a positive constant 
tq such that when t < tq, the time- discrete system i3.1\) - (3m ) admits a unique 
solution (i-*", {7",C"), n = 1, • • • , A^, which satisfies 

WP'^Wh^ + WWm + \\n\w^. + ||AC"||l^ + ||vr|Uco + ( ^ r||VA^7"+ii3/2 < C, 

(3.5) 

and 

max ||e"||Ls + max ||e"||Ls + max ||e"||L'' + f T||Ve"|| L I < Cr. (3.6) 

l<n<Af ^ l<n<7V l<n<Af ^ \^ M c n/, / 

Proof. It suffices to establish the estimates presented in (|3.5p . With such estimates, 
existence and uniqueness of solution follow a routine way. 

We observe that ep, e" and e" satisfy the following equations 



V. ( 4^Ve"+i 1=V- 



k{x) fc(x)\^ 



(3.7) 



$Aee+^ - V • (Z)([/"+i)Ve:?+i) + [/" • Ve^+i 

= V • f (Z)(C/"+^) - D(u"+i))Vc"+^) - ([/" - u") • Vc"+^ - e:?+^/ + 



(3.9) 

for X G i7 and t G [0,T], with the initial and boundary conditions 

^^"^^ -Ve"+i-n = 0, Z)([/"+i)Ve^+i • n = 0, for x G 517, tG[0,r], 



e°(x) =0, for X G , 

(3.10) 



6 



where 

is the truncation error due to the discretization in the time direction. From the 
regularity assumption for c in (j2.ip . we see that 

II^^+^IIl^ < E < ct\ (3.11) 

n=0 

To prove the error estimate p.6p . we multiply p.7p by e""^^ to get 



l|Ve^+i||i. < 



A;(2;) k{x) \ ^^n+i 



which together with (j3.8p implies that 



<C||e^||i2||Vp"+iL- <C||e^||i2 

L2 

(3.12) 



||er^||^2 <C||Ve^+i^.+C||erii.||Vp"+ii^ <C||e^||^.. (3.13) 
Then we multiply (I3JD by ej?+^ to get 

iA($||erii.) + ||V^(C/"+i)Verii. 

< C||e"+i 11^4 Ik^ - + C\\elU. II 1 1|^. || Vc"+i 

< C\\e:^'\\% + C7|||eri||Hi(||erii2 + ||e^||iO + C||f"+i||^_i 
<e||Veri||i,+Q||eri||i.+C7||e^||i,+C|r+i||^_„ (3.14) 

where we have used the inequalities 

|(?7-+i • Ve^i, 6^1)1 = |(V |eri|2)| < ||e-+ii4||/ - g^H^^, 

and 



e 



"+^iii.<^iiveri|ii.+aiieri|ii. 



The square root of D{U^^^) exists because D(C/"^^) is a symmetric and positive 
definite matrix. With Gronwall's inequality and (j3.1ip . (|3.14p reduces to 



N-l 



n=0 

which further produces 

llAe^ii^ <C. 
From ([3T^ - (IXT3] ). we derive that 

l|eML2 + ||Ve^+i||i. <C7r. (3.16) 
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To prove (j3.5p . we rewrite ()3.9p as 



^Dte 



n+l 



V- (L>(u"+i)Ve"+i' 



-V • ((i?(u"+i) - i?(C/"+i))Ve^+i 
+ V • ({D{U"+^) - D(u"+i))Vc"+i 



u 



Ve' 



n+l 



u" • Ve 



n+l 



- u") • Vc"+i - e^+ig^ + 



(3.17) 

Since u*^ G C"(r2), using the W^'^ estimates of elliptic equations for Cc 0, we get 



llVe^iie < C||Aer^||i2 +C7||eriL^||Verii6 + C||u"||i^||Verii. 

+ C||eri||L^||Vc"+i||^6+C||(?^||^6||eri||^3+C||f"+i|U2 



which further reduces to 

||Ve-+i^6 < Co||eri||Loo(Co + HVe^iUz^s) + Co (3.18) 

for some positive constant Co independent of n, r, /i. 

On the other hand, we rewrite the equation ()3.7p into the following form: 



Since we have 



V 



V 



fc(x) 



/i;(x) k{x) 



/i(C«) ;u(c") 



n+l 



(3.19) 



< 



L6 



V 



k{x) 



L6 



n+l I 



<(c + c||vc"||i6)||ve:; 



n+l I 



and 



V 



k{x) k{x) 



< 



V 



/ k{x) k{x) 



L6 



+ 



L6 



A;(x) /c(x) 



V-Vp 



n+l 



L6 



n+l I 



< ||Ve^|U6||Vp"+iL^ + ||e^||i^||p''--|lH'2.«, 
by the VF^'^ estimates of elliptic equations [7], we derive from (I3.19|] that 

\\e;^^w^,e < (C + C||VC"||i6)||Ve^+iLoo + C||Ve^||i6 + C\\e^\\Lo^ 

< (C + C||Ve^||,.6)(6||e^+iH/2,6 + CJepll,.^) 
+ (C + 6)||Ve^|Ua + (C + C,)||e?||,.2 

< (C + C||Ve^||i6)e||e^+i||H^2,6 + C, + CjVe^||i6 (3.20) 



From (13.81) we observe that 



K+'\\w^,.<iC + C\\e^Jw^.e){\\Ve. 



n+l I 



+ l|Vp' 



n+l I 
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<{C + C||Ve^||i6)e||e^+^||i^2,6 + + CjVe^llie, (3.21) 
and by the Sobolev interpolation inequality, we have 

With the estimates (j3.18p , (j3.20p , (j3.2ip and (j3.22p , we now apply mathematical 
induction to prove 

||VeJ?||i6 < 4Co (3.23) 

for n = 1, 2, • • • ,N . Clearly, the above inequality holds when n = 0. With the above 
inequalities, by choosing a proper small e = e(Co), we derive from (j3.20p - (|3.22p that 

\K-^^\\l^<Cit^'\ (3.24) 

where Ci may depend on Cq. Hence, there exists a positive constant tq = tq{Ci) 
such that when r < tq, we have 

Co\K+^\\l- <l/2. 
Substituting the above inequality into (|3.18p gives 

||Ve^+i||i6 <4Co. 

By mathematical induction, (j3.23p holds for 1 < n < A^. From (j3.20p and (j3.2ip . 
we also get 

max 1 1 ep"*""^ 11^4/2,6 + max He""''"'" Hi^i.e < C. (3.25) 

1<71<A'' l<n<N 

By a similar approach, we have further 

||VAf/"+'|lL3/2 < ||VAe::+ii3/2 + II VAU"+^ 11^3/2 

< II Ve;:+i 11^3/2 +Cr-i II Ve;: 11^3/2 + || VAu'^+iL3/2 

< Cr-i||Ve^||i3/2 +Cr-i||Veri||^3/2 + || VAu'^+iL3/2 , 

and therefore 

{^r|VAt/"+ii3/2) <C. (3.26) 

With the above estimates, by applying the estimates of elliptic equations 
[H], we can derive that 

\K^'\\H^ + \K'^'\\m<C. (3.27) 
For the regularity of C""*"^, we rewrite ()3.3p as 

-V • (D ([/"+! )VC"+i) + C/"+i • VC"+i = -$AC"+^ + - r+^q^. (3.28) 
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With llDtC""*"^ ll^^a < C, we can apply the estimates for the ehiptic equation [16] 
to obtain 

||C"+i||^2 < C. (3.29) 

Note that ^ for some a > 0. With the Holder regularity of C", applying 
the VK^'* estimates to ()3.7p and using (jS.Sp . it is not difficult to see that 

lle^^L^ <C||e^||L.. (3.30) 

Multiplying (|3.9p by |ej?"'"^|'^~^e"'*'^ and using (j3.30p . one can derive that 

/ s-i$A|e^+Vdx + (s-l) / \e^+^\'-^D{U''+^)Ve'^+^ ■Ve]^+^ dx 
Jn Jn 

< f f s-i(g^-/)|e^+Vdx+||Vc"+^||Loo||[/'^-u"||L.||e^+^ii^-i 



Jn Jn 

3-2 

+ (s-l)||L»(C/"+i)-D(u"+i)||L.||Vc"+iL.o||e^+i||g^|||e^+V"'Ve^+^||i2 

7n+l||s ^^—InU^+l 11^ _L lU"-!!* ^_L^('£,_^^ / U"+l|''~2|Y7^n+l|2 , 



< (C + ||£:"+l^H^-,.)e-^(||er^||i. + lle^lliO + e(. - 1) / \err~'\^errdx. 



Choosing a proper small e and using Gronwall's inequality lead to 

K+'\\l'<Ct, 
where we have noted the fact that 



l^^+^L^ <C, ^II^^^Tw^-i.^ < Ct\ (3.31) 



n=0 

It follows that 

llAe^+^lli^ <C. (3.32) 

With the above estimate, a refinement of the above arguments shows that 

||ej?+iiy2,. < C7 (3.33) 

and by the Sobolev embedding theorem, ||Ve""'"-'^||Loo < C||e"+^||^^2,s < C. 

Finally, by applying the VF^'* estimates to the elliptic equation (13. 7p and using 
p.Sp . we obtain 

\\'^e;+^\\Ls + \K+^\\Ls <Ct. (3.34) 
The proof of Theorem 13.11 is complete. I 
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4 Error estimates of the fully-discrete system 

To provide optimal error estimates for the fully discrete scheme (|2.6p - (j2.8p . we define 
three projections below. 

Let Hh ■ Lp'iyt) — )• be a projection defined in [34J by 

For any fixed integer n > 0, let 11^ : H^{Q) — )• V/^ be a projection defined by the 
following elliptic problem, 

(z)(?7")V(r; - U^v), = 0, for all £ Vh, v£ H\VL) (4.1) 

and = njj. Moreover, let '■ H — )• be a projection defined in [33] by 

{y -{w- Qhw) , iph^ = 0, for all ^ph £ Sh, w£ H. (4.2) 

By the theory of Galerkin and mixed finite element methods for linear elliptic 
problems [30l[34], with the regularity U"' £ H'^{Q), we have 

||t> - Ilhv\\L2 + h\\v - Uhv\\jji < Ch'^\v\jj2, for all v £ H'^{^1), 

\\v - HI+^v\\lp + h\\v - W^+^vWwhP < Ch'^\v\w2.P, for all v £ W^'P{Q), 

\\w — Qhw\\L'2 + ^11^ ~ Qhwlln'^ ^ Ch'^\w\H2, for all w £ H, (4.3) 

Previous works on Galerkin (or mixed) FEM for the nonlinear parabolic system 
(|l.ip - (|1.3p required the estimate 

\\dt{c^+' -III+^c^+')\\l2 <Ch^ (4.4) 

for an elliptic projection operator H^^^ defined by 

(^D{u)V{v - Uhv), VcPh) = 0, for all 0/, £ Vh, v£ H^Q). (4.5) 

The inequality (j4.4p was proved by Wheeler [39] based on the regularity assumption 
||V£'(u)t ll^oo < C. In the following lemma, we will prove an analogues inequality: 

.N-l . 1/2 

( ^ r|| A(C"+^ - n;j+iC"+i)||^_i J < Ch^, (4.6) 

based on weaker regularity of C/" proved in the last section. The above inequality 
is necessary to obtain optimal error estimates. 

Lemma 4.1 With the regularity assumption ^2.1]) . the regularity o/ (C"""*"^, C/""'"^) 
given in 113. 5|) . and the error estimates given in i3. 6\) . the estimate (^7^ holds. 
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Proof We only prove the lemma for the 3D problem (with s > 3 in Theorem 13. ip . 
The 2D problem can be handled similarly. Note that 

(i?(t/"+i)V(C"+i - n^+iC"+^), V(t>h) = 0, (4.7) 
Z)(?7")V(C"+i - n^C"+i), V^h) = 0. (4.8) 



The difference of the above two equations gives 

_L>(c/"+i)v(n;jc"+i - n^+ic"+^), V(t>h^ 
+ ((L>(c/"+^) - z)(?7"))v(c"+i - n;jc"+i), V(i)h) = 0. 

By the W^'^ estimates of elliptic projections [30], we have 

||v(n;jc"+i - n^+^r+i) 11^6/5 < C7||(L>(c/"+i) - d([/"))v(c"+i - n;:c"+i)||^6/5 

< C||L'(C/"+i) -Z?(?7") 11^2 II V(C"+i -n;jC"+i) 11^3 < C||L»(C/"+i) -!)([/") ||i2/i. 
For any ip £ H^{Q) with J^ipdx = 0, let ^/^ be the solution of the equation 

-V • (Z)(?7"+1)V'0) = if 

with the boundary condition Vip • n = on dQ. Since C/""'"^ is uniformly bounded 
in if^(r2), i.e. HC/""*"-^ ||j|^2 < C, it is easy to check that 

UWh^ < C\\(p\\hi- 

By the boundary condition t/"^^ • n = [/" • n = and the expression of the function 
D{u), we see that D{U"+'^)Vi) ■ n = £>([/") VV' • n = on 917. We see that for 
n > 0, 

(n;jc"+i - n;;+^c"+\ (/^) 

'z)([/"+i)v(n;jc"+i -n;j+ic"+i), v^/;) 

D([/"+i)v(n;jc"+i - n;j+ic"+^), v(V' - p^V^)) 

- ((L>(c/"+i) - D([/"))v(C"+i - n;jc"+i), v(p?,v - V')) 

- ((Li(c/"+i) - D([/"))v(c"+i - n;jc"+i), V?/;) 

< C7||Z)(t/"+i) - I)(C/")||^2/i||V(V^ - P;,V')IIl6 

+ iir+1 - n;jc"+i||^3||v • ((z)([/"+i) - i)(c/"))vv^) ||^^^^ 

<C7||Z)([/'^+1)-I)(C/")||^2/i2||V;||h.2,6 
+ C/i2|r+ivF2,3(||V(I)(C/"+i)-I)(C/"))||^3/2||VV'||L- + ||f/"+^-?7-||^2||V^ 

< C (||AC/"+'||l2 + C||VAf/"+'|lL3/2 + ||VC/"+iAf/"+'|lL3/2) t/i2||V;||h3 
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< C {\\DtU^+^\\L2 + C7||VAf/"+'|li3/2 + ||VC/"+iAf/"+'||L3/2))r/i2||(^||^i. 
Therefore, 



-1 



< C(||Af/"+'||L2 +C||VAC/"+'|Il3/2 + ||VC/"+iAf/"+'|lL3/2))r/i2. 

By (I3.6p . we have 

||AC/"+iL2 < llAe^'llL^ + ||Au"+i,.2 < Ct + IIAu'^+'IIl^, 
||VC/"+1A^/"+^|Il3/2 < ||V?7"+i||i6||A^7"+iL2 < Cr+ ||Au"+^||i2. 

With the regularity assumption ()2.ip on u and the estimate p.6p . we derive that 



Since 



||A(C"+^-n^+'c"+i)l|jf-i 

< II AC"+^ - n^AC"+ij^-i + r-i||n;j+ic"+^ - n^c^+^n^-i 

< c|| AC"+ii^i/i2 + r-i||n^+ic'"+^ - n^c^+i^-i 

< c||Ae:?+i/^i/i2 + c||Ac"+^lbi/i^ + T-i||n';^+ic"+i - liffi'^+^WH-i, 

with the estimate p.6p and the regularity assumption (|2.ip . ()4.6p follows immedi- 
ately. The proof of Lemma l4.1l is complete. I 

Theorem 4.1 Suppose that the initial-boundary value problem l^l.l\) - (L^ has a 

unique solution {p, u, c) which satisfies h2. Then there exist positive constants Hq 
and To such that when h < Hq and r < tq, the finite element system h2. 6\) -i2.'/) 
admits a unique solution {PJ^, UJ^,Cf^), n = 1, • • • , A^, which satisfies 

max IIP." -P"|U2 + max ||C/? - C/"||r2 + max HC - C"|| r2 < C/i^. 
l<ri<Ar" " "^0 l<n<N" " " l<n<N " " " ~ 

where (-P", n = 1, • • • , A^, are i/ie solution of the time-discrete system 13. 



Proof. The solution to the time-discrete system (|3.ip - (|3.4p satisfies 

'-U^^\vh]=-(P^+\V-Vh), (4.9) 



V-U^+\^h) = f/-g^, V^O, (4.10) 
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+ (U^ ■ VC"+\ = (cq' - C^+'q'', ^f,) . (4.11 



for any G Hh, '■Ph G Sh and ^/i G V/j. The above equations with the finite element 
system ([illl-dliHl) imply that the error functions - n/^P^+i, [7^+^ - [7"+i, 



V . (c/;;^+i - [/«+!), c^^) = 0, (4.13) 
$A(c^' -c"+i), 0;,) + (z)([/^"+i)v(c;^+i -n-+ir+i), v</.;,) 

(q+i _c"+i)g^, 0,) + ((i)(c/"+i) -i)(c/,"+i))vn;:+ic"+\ v^,). (4.i4) 

First we present a uniform estimate for the finite element solution. Let g be the 
solution of the elliptic problem 

with the boundary condition -^0r^^g • n = on dVl. 

k{x) 



Taking Vh = Qh [ u(c" ) ^9 1 P-Sp and by noting the definition of the projection 



Qh, we get (with Vh ■ n = -^^Vg ■ n = on d^l [34]) 



<c\\uj:^'\U\vHy 
<c\\ur'\uw-<i''\\L^- 

and taking Vh = UJ^^^ in ([221) and = ^'/T"^^ in gives 

= -(^^^ v • f/r^) = / - / 

Combining the last two inequalities, we obtain 

<C||<7^-g^||^,. (4.15) 

Secondly, we present an upper bound for the error function \\UJ^^^ — U"'~^^\\i2 in 
terms of HC^"*"^ — C""'"^||/^2. By the definition of the projection operator Qh in (j4.2p . 

(X/ . ([/^"+i - QhU^+'), ^h) = 0, for all G S^, 
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which together with (|i33]) shows that V • (C/^+^ 
Vh = - QhC/"+^ in we get 



QhU 



in ri. Taking 



k{x) 



QhU'^^')+^{QHU^+'-U^+') 



+ 



\ k{x) k{x) 



jjn+l^ jjn+l 



QhU 



n+1 



0, 



u 



n+1 1 



L2- 



which further imphes that 

Wu-^+^-QhU^+^W^.KCWQhU^^^ 
With ()4.3p . the above inequahty reduces to 

C/"+li2 < Ch^ + C\\Cl-C'^\\^^. (4.16) 
/^(C^+i - Iil+^C^+^)dx (so that 



I 



Thirdly, we take (j)h 
(ph dx = 0) in (|4.14p and get 



jjn+l^n+l 



2 



<6||V(C 

+ C\\q' 



$(c;^+i -n;j+ic"+i) 

n+1 ■Qn+l^n+lM|2 



L2 



nr'c"^')lli2 + ||A(c 



Il3|IC 



+ A/D([/r^)V(C^ 

nn+l 



n+1 j-|-n+l^n+l 



L2 



->n+l 



n+1 1 



, -n;j+ir+i)iu.||cr^-c"' 11^2 



L- IIL2 
n+1 1 



+ c||c/"IIlo.||v(c;^+i 
+ c||c;^+i -n;^+ic"+ii6(||?7^"+i - c/"+^||i2||v(n;:+ic"+i)ll^3 
+ - c/"+^||l3||v(c^i - ui+'c-+')\y) 
+ C|k^||i3(|ic+i-n;^+ic"+i||i3 + ||c;; 



nr'C"+^llL6||c 



+ c'll^n^'^ c ■ iiLo^ii V 



< e||V(C' 

+ c||c;^+i-n 



nn+l(.n+l^ 



- u 

2 



-in+l T-rn+l 

- '-'^h IIL3 -r ||L 

jn+l^n+1 
,n+l _ TT"+lcn+l)||2^ + || A(C"+^ 

- f/"+i^2||V(C"+l -ff;^+lC"+l)||2 2 

where we have used (j4.3p . (j4.16p and the following estimate: 



n+1 



nr^c-+i||^2) 



n+1 1 



|L2 



+ C/l-l/2||f7n+l 



(4.17) 



jjn+l^n+l 



+ 



f^n+l^^n+l 



^n+1^ ^n+l -Qn+l^n+1 

C"+i), V(C, 



->n+l 



jjn+l^n+1- 



n+1 1 



IL2 



->n+l 



nn+l^n+l^ii^^liq+l _^n+l|| 



L2 



Prom (|4.16p we observe that (j4.17p reduces to 

2 



jjn+l^n+l) 



+ 



L2 



Z)(t/^'+i)V(C^ 



L2 
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+ c\\ci^^ - n^+^r+^ii^^ + c\\ci - TilC'Wl^ + II A(C"+i - n;;+ic"+i)||5^_i + ch^. 

By applying Gronwall's inequality with induction, we deduce that there exist /iq > 
such that for h < ho, 



■n+l 



By ()4.16p . we get further 



(4.18) 

(4.19) 
(4.20) 

Finally, we estimate the error \\Ph — P"~'~^||^2. We redefine g to be the solution 
to the equation 



ll^n+l < Ch^, 



with the boundary condition -^feryV^ • n = on dQ. Easy to check that 



Let 



k{x) 



vh = Q 



rV5 



Then 

{ip,V • Vh) 
and from (j4.12p we obtain 



(v9,p^"+i-n;,p"+i), ^GSh 



,n+l||2 
IIl2 



A;(x 



/c(x) 



A;(x) 



-V5 



+ 



\ A:(x) /c(x 
k{x) ^ 



-V5 



k{x) 



+ 



MQ) 

A;(x) k{x) 



k{x) 



where we have used ()4.13p . By (I4.16p . we have further 



HI 
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+Ch^p ■ - f/"+i)|L.||5||H. + C\n - C"||L2||C/,"+i||^3||V<7||i6. 

Since 

I|V5||l6 < ciblln^ < c||p,"+i -n;,p"+i||i2, 

and 

< \\g\\H2 + \\VC^+^\\lA\9\\h^ < Cll^r' - n/.P"+^||L2 , 

//I 

it follows immediately that 

The proof of Theorem 14.11 is complete. I 

Combining Theorems 13.11 and Theorem 14.11 we complete the proof of Theorem 

EH ■ 

5 Conclusions 

We have studied error analysis for a nonlinear and strongly coupled parabolic system 
from incompressible miscible flow in porous media with a commonly-used Galerkin- 
mixed FEM and linearized semi-implicit Euler scheme. Optimal L? error estimates 
were obtained almost without any time-step (convergence) condition, while all pre- 
vious works require certain restriction for time-step size. We believe that the idea 
of the error splitting coupled with the regularity analysis of the time-discrete PDEs 
can be applied to many other nonlinear parabolic PDEs and time discretizations to 
obtain optimal error estimates unconditionally. 
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